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(N ■ Abstract 

We study formulation and probabilistic interpretation of a simple general-relativistic hamiltonian 
quantum system. The system has no unitary evolution in background time. The quantum theory yields 
. transition probabilities between measurable quantities (partial observables). These converge to the clas- 

sical predictions in the h — > limit. Our main tool is the kernel of the projector on the solutions of 
Wheeler-deWitt equation, which we analyze in detail. It is a real quantity which can be seen as a prop- 
agator that propagates "forward" as well as "backward" in a local parameter time. Individual quantum 
states, on the other hand, may contain only "forward propagating" components. The analysis sheds some 
light on the interpretation of background independent transition amplitudes in quantum gravity. 

o ' 

rn ! 1 Introduction 
O 

Conventional mechanics describes evolution in "background time": the independent time variable t is 
assumed to represent a measurable but non-dynamical physical quantity. This is true in the classical as well 
as in the quantum theory, where evolution in background time is given by unitary transformations. With 
general relativity, however, we have learned that there is no non-dynamical measurable time in nature: there 
is no background spacetime, and no background time in particular. Therefore, at the fundamental level 
physics cannot describe evolution in time. It can only describe relations, or correlations, between measurable 
quantities. It is therefore necessary to extend the formalism of classical and quantum mechanics to such a 
background independent context. 

A strategy in this direction has been explored in a number of recent works ^ |2J E] 1 Here we study 
a simple system in order to illustrate and test ideas and techniques discussed in these works. We do not 
restate the general theory here, for which we refer the reader to the references above, but this paper is 
nevertheless self-contained. The system we consider has one degree of freedom. It has a well defined 
classical dynamics. It could describe a simple cosmological model - or a simple mechanical system studied 
without using an external clock. The system has two variables, which we call a and b. Locally, we can view 
b as the independent variable and a as the dependent one. The model describes then the evolution a(b) 
of a as a function of b. In other words, locally we can view b as the time and a as the physical degree of 
freedom. But globally this picture breaks down: the dynamics predicts relations betweent the two variables 
a and b which are on an equal footing. The quantum evolution in b fails to be unitary. This is the difficulty 
that generates the much debated "problem of time" in quantum general relativity [7] . 

The system we consider has been investigated in (51 El- The difficulty is not to define the quantum 
states of the theory - a rather easy task using standard hamiltonian methods. The difficulty is to define 
a consistent probabilistic interpretation of the formalism. Due to the lack of unitary time evolution, the 
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conventional probabilistic interpretation of the wave function at fixed time does not apply. A good discussion 
of the difficulties of the "naive" interpretation of the wave function as a probability density, and of tentative 
strategies to solve the problem can be found in In one of us has studied the physical interpretation 
of this system using operators well defined on the physical Hilbert space. These describe "complete" 
observables, (or "perennials") and include, in particular the "evolving constants of the motion" 0IE]- This 
way of interpreting general relativistic quantum systems is correct in principle, but may be cumbersome in 
practice. 

In PP, it was pointed out that the conventional notion of "observable" may not be the most suitable one 
in a background independent context. The weaker notion of "partial observable" may be more useful in such 
a context. A partial observable is any quantity that can be measured - even if it cannot be predicted by 
the theory. Using this notion, general relativistic systems admit a simple reading: they predict correlations 
between partial observables. This is the point of view developed in 00 In this philosophy the key object 
that yields the interpretation of the quantum formalism is the "propagator" K(a,b,a' ,b'). This is defined 
as the kernel of the projector from the kinematical Hilbert space /C (the space representing all possible 
outcomes of measurements of partial observables) to the physical Hilbert space H. (the space representing 
the sole outcomes of measurements of partial observables that are allowed by dynamics). The propagator 
plays a role similar to the propagator of the Schrodinger equation in nonrelativistic quantum mechanics, 
which, as emphasized by Feynman, expresses the quantum dynamics. The propagator K(a,b,a' ,b') can be 
given a direct physical interpretation This point of view might prove useful for the interpretation of 
the formalisms of non perturbative quantum gravity - in particular in loop quantum gravity |Hj and its 
sum over path version, the spinfoam formalism ^U], which gives us the analog of K(a, b, a' , b') in quantum 
general relativity 

Here we study the propagator K(a, b, a', b') of our simple model and we discuss physical interpretation 
and classical limit. We show that in a regime in which b can be taken as an independent time parameter 
the probabilistic interpretation proposed in [2] reduces to the conventional interpretation of the modulus 
square of the wave function as a spacial probability density at fixed time. We illustrate the semiclassical 
limit and some semiclassical states. The propagator turns out to be real. It can be seen as propagating 
"forward" as well as "backward" , in a local parameter time, in spite of the fact that the difference between 
forward and backward propagation is not manifest in the classical theory. The quantum states, on the other 
hand, split among forward and backward propagating ones. This might be the same situation that we find 
in quantum general relativity. 

In the section 2 we discuss the classical theory. We formulate the theory in a fully covariant language, 
that illustrate the structure of general-relativistic hamiltonian mechanics. We study the quantum theory in 
Section 3. In the appendices, we recover an approximate Schrodinger equation, we recall the probabilistic 
interpretation based on the flux of the conserved current, and we illustrate a derivation of the propagator 
using functional integral techniques. A complete formulation and a general discussion of the classical and 
quantum formalism that we utilize here can be found in chapters 3 and 5 of reference 

2 Classical theory 

2.1 Definition of the model 

A simple harmonic oscillator can be viewed as system with two partial observables, q and t. A motion of the 
system defines a relation between q and t. A given motion is characterized by the two constants A £ [0, oo] 
and (p £ [0, 2ir], and is given by the equation 



where u is a constant characterizing the system. This motion is a sinusoid of amplitude A and phase <j> 
in the (q, t) plane, the extended configurations space C. This dynamics is the classical limit of a quantum 
dynamics. 




(1) 
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Here, we consider a system with two partial observables a and b. A given motion of the system is 
characterized by two constants A £ [0, y/2M] and (j) G [0, 2tt], and is given by the equation 

f(a, b) = (£) 2 + 2 - 2 cos ~ = sin 2 0, (2) 

where M is a constant characterizing the system and B 2 = 2M — A 2 . This motion is an ellipse of radii A 
and B and inclination (j) in the (a, b) plane, the extended configuration space C. We want to understand 
if this system can be seen as the classical limit of a quantum system, with a well defined probabilistic 
interpretation. 



2.2 Classical dynamics 

A relativistic hamiltonian formulation of the system above can be given on the extended configuration space 
C = R 2 coordinatized by the partial observables (a,b). The dynamics © is governed by the relativistic 
hamiltonian 

H(a,b,p a ,p b ) = l(a 2 + b 2 +p 2 a + pl)-M (3) 

defined on the cotangent space T*C = R 4 coordinatized by (a, b,p a ,p b ), where p a and p b are the momenta 
conjugate to a and b. In general, in a general-relativistic hamiltonian system the dynamics is defined by 
the surface H(q a ,p a ) = 0, in T*C (see chapter 3 of 00 The hamiltonian (J3J) yields the Hamilton-Jacobi 
equation 

1 ' + — ^t^- + a 2 + b 2 - 2M = 0. (4) 



V da J \ db 

Given a one parameter family S(a, b; A) of solutions of this equation, the physical motions are determined 
-given two constants A and p^— by the equation 

dS(a,b;A) 

OA ~ PA = - (5) 
A family of solutions of the Hamilton-Jacobi equation Q is given by 

A 2 ( a \ b f-= B 



S(a, b; A) = -y/A 2 - a 2 + — arctan - + -VB 2 - b 2 + — arctan , . (6) 

V ' 2 V ^2 \^A 2 -a 2 J 2 V ^2 VV-B 2 - b 2 J V ' 

inserting this in © gives @, where pa = Acft, with simple algebra. 



2.3 Action 

Consider a path A : [0,1] — ► T*C in the space of the coordinates a, b and their momenta p a ,Pb, and let 
A : [0, 1] — * C be its projection on the extended configuration space C (just the coordinates a, b). The action 
S [A] is defined by 

S[X) = J ( Pa da + Pb db). (7) 

The physical motions A are the ones such that A extremizes this action in the class of paths that satisfy 

H{a,b, Pa , Pb ) = (8) 

and such that A has fixed extrema A(0) = (a, b) and A(l) = (a',b'). These motions follow the ellipse (|2|. 
This is a very general structure, common to all non relativistic as well as general relativistic systems 0. 
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2.4 Geometry 

Since the space T*C is a cotangent space, it is naturally equipped with the Poincare one-form 

9 = p a da + p b db (9) 

and the symplectic two- form uo = dO. Notice that the action is simply the line integral along the path 
of the Poincare one-form 

S{\] = ( 9. (10) 



Equation (JHJ) defines a surface £ in the space T*C. In the model we are considering this surface is a three- 
sphere of radius y/2M. The restriction of uj to £ is degenerate: it has a null direction X. That is, there is 
a vector field XonE that satisfies uj(X) = 0. It can be easily verified that this is 

9 d d , d 

oa do opa opb 

The integral lines of X on £ define the motions. It is easy to see that these integral lines are circles. In 
fact, they define a standard Hopf fibration of the three-sphere 5*3 ~ S2 x Si. The motions are precisely the 
ones given by @, and the corresponding momenta are 

p a = \J A 1 - a 2 , 

Pb = \/B 2 - b 2 . (12) 

The quotient of S by the motions S2 ~ S3 /Si is the physical phase space of the system, coordinatized by 
the two constants A and 4> 

2A 2 = 2M + p 2 a - p\ + a 2 - b 2 (13) 
. , p a b-p b a 

tan0 = ■ — -, (14) 

PaPb + ab 

which are constant along the motions. 

2.5 Parametrization and the associate nonrelativistic system 

The dynamics of the system becomes easier to deal with if we introduce a non-physical parameter r along 
the physical motions (J2J). We can parametrize the ellipse (J2J as follows 

a(r) = A sin(r + tfi), 

b{r) = B sin(r). (15) 



The momenta are then given by 

p a {r) = \/ A 2 — a 2 = A cos(r + 



da(r) 



dr 



db{r) 



Pb (r) = VB 2 - b 2 = B cos(r) = =^-. (16) 

or 

The parameter r is not connected with observability; it is only introduced for convenience. In particular 
the physical predictions of our system regard the relation between a and b, not the dependence of a and 
b on r. The parameter r corresponds to the time coordinate of general relativity: physically meaningful 
quantities are independent from the time coordinate. 
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It is useful to consider also a distinct dynamical system with two degrees of freedom a(r) and 6(t) 
evolving according to equations (|15|) in a time parameter r. We call this system the associate nonrelativistic 
system. It is important not to confuse the two systems. The associate system has two degrees of freedom, 
while our original general-relativistic system has one degree of freedom. The associate system is governed 
by the function H given in Q, seen now as a conventional non-general-relativistic hamiltonian. Indeed, 
it is immediate to see that equations (|15I16|) are the solutions of the Hamilton equations of H. This is 
the hamiltonian of two harmonic oscillators with unit mass and unit angular frequency, minus a constant 
energy M? 

The original system can be derived from the associate nonrelativistic systems in two steps. First, by 
restricting the motions to the ones with total energy equal to zero (namely the ones in which the energy 
E = E a + E b of the two oscillators is E = M). Second, by gauging away the time evolution in r: that is, by 
restricting the physical observables to the r independent ones, namely to the sole relations between a and 
b. This reduces by one the number of degrees freedom. Intuitively, we "throw away the clock that measures 



2.6 The Hamilton function 

As beautifully emphasized by Hamilton ^2], the solution of any dynamical system is entirely coded in its 
Hamilton function. This is true for background independent systems as well. Furthermore, the Hamilton 
function plays a key role in relating the classical system to the quantum system, and we shall use it 
extensively. It is essentially the classical limit of the quantum propagator. 

The Hamilton function S(a, b, a' , b') is a function onCxC that satisfies the Hamilton-Jacobi equation in 
both sets of variables. It is defined as the value of the action of the physical motion that goes from (a' , b') 
to (a, b). If S(a,b,a' ,b') is known, all physical motions can be obtained just taking derivatives, as follows. 
Defining 

, , , dS(a,b,a',b') , , , dS(a,b,a' ,b') 
p a {a,b,a,b) = — , p b (a,b,a,b) = — , (17) 

the algebraic equations 

p' a (a, b, a', b') = p' a , p' b (a, b, a', b') = p' b (18) 

are equivalent to © and give the physical motions. Here a motion is characterized by the "initial data" 
(a', b' \p' a ,p' b ) instead than by A and (ft. 

Let us study the Hamilton functions of our system. Given two points (a',b') and (a, b) in C, we first 
need to find the physical motion that goes from one to the other. That is, we need to find the value of the 
two constants A and (ft of this motion. To this aim, it is convenient to parametrize the motion as in ()15l) 
and search the value of the three constants A, (ft and r such that 

j a' = o(r') = Asm(T ! + (ft), ja = o(r' + r) = A sin(r' + r + (ft), . , 

) b' = b(r') = B sin(r'); \ b = b(r' + r) = BsimV + r). 1 j 



With some algebra we find 



and 



■ 2 1 a' 2 — 2aa' cost 
sin 2 r 



A 2 = — — (20) 



(a 2 + b 2 + a' 2 + b' 2 ) - 2{aa' + W) cos r 
M — — 2 (21) 
sin t 



2 In turn, the associate system can be cast in covariant form as well. Its extended phase space has coordinates (a, b, r) and 
its relativistic hamiltonian is H associate — Pt + H. 
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The last equation can solved for r, giving 

, , , lA aa' + W ± J(aa' + bb'Y + M(M - a 2 - b 2 - a' 2 - W) , , 

t(o, 6, a', 6') = arccos Z± - — —± - (22) 

Inserting this value in (j2U|) gives ^4 = A(a, 6, a', b'). 

The value of the Hamilton function is obtained by performing the integration that defines the action 
(|7|) of the motion that goes from (a', b') to (a, 6). The result is 

S(a, b, a', if) = S(a, b, a', b'; A{a, b, a', &')) (23) 

where 

S(a, b, a', b'; A) = S(a, b, A) - S(a, b', A). (24) 

Here S(a,b,A) is given in © and A(a,b,a' ,b') is the value of A of the ellipse (J2J that crosses (a,b) and 
(a',b'), determined by the equations (|2*U)) and (|2*2*|) . 
It easy to show that 



dS(a,b,a',b'; A) 



OA 



= (25) 

A=A(a,b,a',b') 



We can also write the hamilton function in another form, that turns out to be useful in the following. Using 
(|15|) and 1)21(1 we can write 

S(a, b, a', b') = S(a, b, a', b'; r(a, b, a', if)) (26) 

where 

\ 1 (a 2 + b 2 + a' 2 + b' 2 )cosT -2(aa' + bb') w 

S 1 a, 6, a', 5'; r = ^ + Mr. (27) 

2 sinr 

This result is not surprising. The Hamilton function of a single harmonic oscillator is 

, v l (a 2 + a /2 )cosr-2aa / 
S(a,a;r) = - : ; (28) 

2 S1I1T 

equation (|27|) shows that the action of the relativistic system is equal to the action of the associate system 
of two harmonic oscillators a(r) and 6(r) plus the constant M term, calculated for the time r determined 
by the requirement that the motion reproduces the constrained one of the covariant system. That is 

S(a, b, a', b') = S{a, a'; r(a, b, a', b 1 )) + S{b, b'; r(a, b, a', b')) + Mr(a, b, a', b'). (29) 

Finally, it is easy to check that 

dS(a,b,a' ,b'; r 



Or 

An equation that plays an important role below. 



= (30) 

T=r(a,b,a' ,b') 
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3 Quantum theory 



The classical system described above can be viewed as the classical limit of a quantum system. The quantum 
system is defined on the kinematical Hilbert space K. = L2[R 2 , da db] formed by the wave functions ip(a,b). 
The partial observables a and b act as (self-adjoint) multiplicative operators. The dynamics is defined by 
the Wheeler-deWitt equation 

( _?l2 £ " **W + ° 2 + 62 " 2M ) ^ 6 ) = °" ( 31 ) 

The operator in parenthesis is the Wheeler-deWitt operator H. The space of the solutions of this equation 
form the Hilbert space 7i of the theory. It is easy to find the general solution of this equation. The equation 
can be rewritten as 

H(a,b,p a , Pb ) VM) = (H a + H b — M) ^(a,b) = 0. (32) 

where H a (and Hb) is the Shrodinger Hamiltonian operator (with unit mass and unit frequency) in the a 
(resp b) variable of a single harmonic oscillator. Since the eigenvalues of this operator are E a = U(n a + 1/2) 
(and Ef, = Ti(rib + 1/2)) with non negative integer n a (and rib), m the basis of the energy eigenstates, \n a ,rib), 
equation (|3l"j) takes the form 

H(a,b,p a ,p b ) \n a ,n b ) = (h(n a + n b + 1) - M) \n a ,n b ) = 0. (33) 

This equation has solution only if 

M = h(n a + n b + 1) (34) 

namely if N = M/h — 1 is a nonnegative integer. We shall assume so from now on. The general non- 
normalized solution of (|3*T|) is then easily 

i^(a,b) = C n a ,n b ^n a (a) ipn b (b) (35) 

n a +n b +l=N 

where ip n is the normalized n-th eigenfunction of the Harmonic oscillator, and the coefficients c na ^ nb satisfy 
J2n a +n b +i=N \ c n a ,n b \ 2 = 1- The state space fC is the Hilbert space of the associate nonrelativistic system. A 
basis in fC is formed by the states \n a ,rib) with n a a-quanta and rib 6-quanta: 

(a,b\n a ,n b ) = ipn a (a) ipn b (b) (36) 

The physical Hilbert space H is the one spanned by the states with N quanta. It is convenient to write 
j = N/2 and to introduce the quantum number m = \{n a — rib) that runs from m = — j to m = j. So that 
H is spanned by the (2j + 1) states 

\m) = \j — m,j + m). (37) 

That is 

(a, b\m) = ipj- m (a) tpj+mip) = ip m (a, b). (38) 

Explicitly, 

ip m (a, b) = = H j+m {a/VK)H^ m {b/y/n) e"T (39) 

y2 i 3h-K (j+m) ! (j-m) ' 

where H n (q) is the n-th Hermite polynomial. (As pointed out by Schwinger we can define angular momen- 
tum operators on the Hilbert space of two oscillators; the total Hamiltonian is proportional to the total 
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angular momentum. TL is then the Hilbert space of the irreducible spin-j component of tC. m is the quantum 
number associated to one of the components of the angular momentum: 

L» = \(vl-pl + a 2 -b 2 ) (40) 

see jH].) Alternatively, we can diagonalize the angular momentum 

L = \{b Pa -ap h ). (41) 

Its eigenstates are easily written in polar coordinates a = r sin <p, b = r cos ip as 

zMr !¥ = e! m '^m'(r), ( 42 ) 

where m! is the eigenvalue of L and ip m i (r) is the (unique) solution of the radial equations with the correct 
energy 

Vw(r) = e~ r2/2 r^'l^-njm'l + iy) (43) 

where n, \m'\ + l,r 2 ) is the confluent hypergeometric function and n = 0,1,..., j, \m'\ = 0,1,..., 2j, 

with the condition 

2j = 2n + \m'\. (44) 

TL is a proper subspace of /C, therefore the scalar product of 1C is well defined in 7i. As well known, 
in several general relativistic systems the space of solutions of the Wheeler-deWitt equation is a space of 
generalized states, namely a subspace of a suitable completion of /C. This is sometimes described as serious 
conceptual difficulty, but it is not: there are many equivalent techniques for deriving a scalar product of 7i 
from the scalar product of K. Here we are not concerned with this technicality. 

3.1 The propagator 

Since TL is a proper subspace of /C, there is an orthogonal projection operator 

P :K^TL. (45) 

The relativistic propagator K(a, b, a' , b') is defined as the integral kernel of P. Easily, this projector is given 
by 

P= ^2 \n a ,n b )(n a ,n b \ (46) 

n a +n b +l=N 

Therefore 

K(a,b,a',b') = ^ (a,b\n a ,n b )(n a ,n b \a' ,b'). (47) 

n a +n b +l=N 

Explicitly, we can write this as 

K(a,b,a',b') = £ y* m (a',b')y m (a,b) (48) 

m=-j 

y> e __k( a 2 +6 2 +a ,2 +6 , 2) H j+m (a/^)H J ^ m (b/Vh)H j+rn (a'/Vh)H j ^ m (b'/Vh) 
4±- 2 2 J'(j + m)!(j - m)!/i7r 1 j 

using the properties of the Hermite polynomials, this can be written for integer j also as 

1 -Hg±^!iV! 4^ {2{aa' + W)/nfi- 2i r2j -2i { a 2 + b 2 + a' 2 + b> 2 \ 
K(a,b,a,b) = -e * g ^— ^ L 4 ^ j (50) 

where is the n-th generalized Laguerre polynomial with parameter m. 
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3.2 Alternative expressions for the propagator 

The expression 1)46(1 of the propagator can be rewritten as follows 

P = X! \n a ,nb)(n a ,n b \ = Y |n , n b )(n a , n b \ 5 na+na+ i- N 

n a +n b +l=N n a ,n b 

= ^ \n ai n b )(n a ,n b \ T dr e^^ +1 ~^ = ^ \n ai n b )(n ai n b 

n a ,n b J0 n a ,n b 

dr e~* H \ 

o 



dr e~% HT 



o 



That is 



2- 



Pip{a,b) = I dr i/)(a,b,r), 
Jo 



(51) 



(52) 



where ip(a,b, r) is the wave function of the associate system that evolves in a time r from the initial state 
tp(a, b, 0) = ip(a, b). The evolution of the associate system is easy to determine. We have 

ip(a, b, t) = e-^ HT i(j(a, b, 0) = J da' J db'K(a, a', r)K(b, b'T)e~l Mr if>(a, b, 0) (53) 

where K(a, a'r) is the well known propagator of a single Harmonic oscillator 



K(a,a',T) = {a\e-* HaT \a') 



1 



y/2mh sinT V 2 sin r 



exp ( f (a 2 + a' 2 ) cos r — 2aa' 

\ 2 sin t V L 



(54) 



and the term e r Mt is present because of the constant term in the hamiltonian. Inserting ()54|) in Q53|) we 
find a simple integral expression for the propagator 



r27r I / i 

K(a,b,a! ,b') = / dr ri _,_ fc - — -exp 



o 2-7ri?isinT \27isinT 



(a 2 + 6 2 + a ,2 + 6 ,2 )cosr-2(a / a + 6 , 6) +^Mr (55) 



that is 



K(a,b,a,b , )= / dr exp I —S (a, b, a', 

Jo 27rih suit \n 

where S(a, b, a', b'\ r) is the hamilton function given by 1(27(1 . 



b'-T 



(56) 



3.3 Properties of the propagator 

Let us summarize some of the properties of the propagator. The propagator 

• satisfies the Wheeler-deWitt equation ((31(1 in both sets of variables. 

• satisfies the composition law 

da' J db' K(a, b, a', b') K(a', b', a", b") = K(a, b, a", b"), 

this follows immediately from the fact that P is a projector; 
propagates a physical state from (a', b') to (a, b) in the sense 

da'db' K(a,b,a' ,b') ip(a',b') =ip(a,b); 



(57) 



(58) 



(59) 



• projects an arbitrary function ip(a, b) (not a solution of to a solution of 

• is gauge invariant, in the sense that it is independent from the parameter r; 

• is real; 

• is symmetric in the exchange (a, b) *-* (a',b'). 

4 Formal relations between quantum and classical theory 

From now we assume that M 3> % and we study the relation between the classical and the quantum theory. 
4.1 The semiclassical limit of the propagator 

We begin by relating the propagator K(a, b, a' , b') and the Hamilton function S(a, b, a' , b'). A semiclassical 
approximation of the propagator can be derived from the expression (|55|). In the limit of small fi we can 
evaluate the integral using a saddle point approximation. The result is 



K{a, b, a', b') = Y - 1 exp ( \(a 2 + b 2 + a' 2 + b' 2 ) cos r< - 2(o'a + 6'6)1 + ^Mr, ) (60) 

. 27ri n sin r,- V 2n sin r,- L J ft < 



where T{ = Tj(a, 6, a', b') are the values of r for which the exponential has an extremum. That is, they are 
determined by 



d 



dr \2?isinr 



(a 2 + b 2 + a' 2 + 6' 2 ) cos r - 2(a'a + 6'6) 



+ -Mr 

h 



(61) 



But notice that the exponent is precisely S(a,b,a' ,b' , r), and equation (|6*T|) is precisely equation (j5U|) that 
determines the classical time r along the classical trajectory from (a', 6') to (a', 6'). Therefore the exponent 
in (|60|) is precisely (jr) the Hamilton function ! 

There is however an essential subtlety, that we have disregarded during the classical discussion. In 
general, equation (|61j) does not has one, but rather two solutions (a smooth periodic function cannot have 
a single extremum). Why two? Because in general there are two paths that connect two points (a',b') and 
(a, b): a shorter one and a longer one obtained subtracting the shorter one from the entire ellipses. Indeed, 
it is clear that if r is a solution, so is 2ir — r. Therefore the sum over Tj contains two terms: t\ = r(a, b, a' , b') 
where r(a,b,a' ,b') is given by (|22j) and T2 = 2tt — r(a,b,a' ,b'). Inserting this in ()62j) we conclude that, to 
the first relevant order in h, 

K(a, b, a', b') = — L— — (Vi^Aa'M _ e+ *S(-MW) (62) 
27rm sm r(a, o, a, o'J V / 

or 



K(a, 6, a', 6') = — - - 7~7/\ sin -5(o, 6, a', 6') . (63) 
7rn sm r(a, o, a , o'J \n / 

4.2 Quasi classical states 

Consider a coherent state of the associate nonrelativistic system. Using standard coherent state technology, 
the state determined by the classical position (a, b) and the classical momenta {p a ,Pb) can be written as 

oo oo ot na (y n ' b ( |q |2 _|_ Iq, |2\ 

\a,b,p a ,Pb) = f i b f exp ^ — 9 — ~ l n «' n fc) ( 64 ) 



n a =u n b 



'o^o Vnjn b \ 
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where a a = a + ip a and a b = b + ipj,. This is a minimum spread wave packet centered in the position (a, b) 
and with minimum spread momentum (p a ,Pb)- Recall that the coherent states of the harmonic oscillator 
have the property 



a,b,p a ,p b ,T) = e k Ht \ a ,b,p a ,p b ) = |a(r), b(r) , p a (r) , p b (r)) 



(65) 



where («(t), 6(r),p a (r),pb(r)) is the classical motion determined by the initial conditions a,b,p a ,p b . That 
is, they follow classical trajectories exactly. 

Let us now project this state on 7i. We obtain 



E 



n a +n b +l=N 



P\a,b,p a ,p b ) -- 
On the other hand, we have also 

ipa,b,p a , Pb (a'i b ') 



y/n a \n b \ 



exp 



IcU 2 + \ a b\ 2 



\n a ,n b ) 



2~ 



dr ipa,b,p a , Pb (a',b',r) 



(66) 



(67) 



where ip a ,b,p a ,p b ( a i t) is the time evolution of the coherent state of the associate system. In the small fi limit 
the state ip(a, 6, r) has essentially support only along a narrow strip around the classical motion determined 
by the initial conditions (a,b,p a ,p b ). Therefore the physical state ipa,b,p a ,p b (a>' \b') will have support on this 
same region, namely on one of the ellipses ((SJ. The figure shows the square of the amplitude of a coherent 
state ip(a, b) computed numerically: it is picked around the classical solution. 




Figure 1: Support of a quasi classical state 



4.3 Forward and backward propagation 

In the classical theory, an ellipse defines a classical motion. We can assume that the versus of the ellipse, 
clockwise or anticlockwise, has no physical significance. In fact the equation © that fixes the motion 
does not contain any reference to the versus of the ellipse. However, in the quantum theory there are 
distinct "clockwise and anticlockwise propagating" states. To see this, consider, instead of the coherent 
state defined in the previous section, the one defined by the classical initial values (a,b,—p a ,—p b )- In 
the associate system, the same ellipses is followed in the opposite direction. After integrating in dr, the 
resulting physical coherent state i/j a ,b,-p a ,-p b has the same support as ip a ,b,p a ,p b , but it is not the same state 
! A moment of reflection shows that it is its complex conjugate 

1pa,b-p a -p b (a',b') = 1pa,b,Pa, Pb ( a '> b ')- ( 68 ) 
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In fact, the physical Hilbert space can be split into the forward and backward propagating subspaces 
according to the sign of the eigenvalues of the operator L. 

A particularly simple semiclassical state, which illustrates well some general features of these states, can 
be obtained as follows. Consider the classical motion formed by a circle, namely the case A = B = y/M. 
The two corresponding semiclassical states are easily obtained. They are the states (|42|) with m' = ±2j, 
that maximize and minimize the angular momentum L. In the large M/Ti limit, the radial function become 
a narrow gaussian around the classical radius y/M, as can be seen from (|43j) . Indeed, when m! = ±2j, 
n = 0, and since 

iFi(-n, \m'\ + l,r 2 ) = 1 (69) 

we have 

^±2j(r,<p) = e -r 2 /2+2ilar e ±i2j>_ (?Q) 

The maximum is in r = ^/2j ~ y/M. 

This observation sheds light on the fact that the propagator is real and contains the two terms in the 
right hand side of (|B*2*|) . The propagator propagates both sectors of the Hilbert space - the one propagating 
clockwise, ip-2j(r,(p), and the one propagating anticlockwise, ip+2j (r, f) ■ Each of the two terms of (jF)2*|) 
propagates one of these sectors. Since the two are simply the complex conjugate of each other, the propagator 
is real. 

On the other hand, there are quantum states that are in only one of the two sectors. They are "clockwise 
propagating" states. Restricted to those states, the propagator can be written as a single exponential of 
the Hamilton function. 

Notice that the two terms of (|62|) can be interpreted in two equivalent manners. Either as related to the 
two branches of the ellipses that connect (a, b) and (a', b'). Or as one term evolving from (a, b) to (a', b') and 
other from (a',b') to (a,b). The math of section l4~T1 can be interpreted either way. Thus, the second term 
of the propagator can be seen as the backward propagating term. We suspect that this term is therefore 
going to be present irrespectively of whether the classical orbits are closed. 

Recall that in quantum general relativity the three-geometry to three-geometry transition amplitude has 
a surprising tendency of turning out to be real. For instance, the Ponzano-Regge transition amplitudes are 
real. Similarly, spin foam model tend to give real transition amplitudes. We suspect that the propagator of 
a covariant theory is naturally defined as a real propagation -as in the example studied here- where there 
is no a priori distinction between forward and backward propagation. On the other hand, a quantum state 
may contain one component only, and therefore be able to "select" the appropriate part of the propagator. 
Form this point of view, the attempt by Oriti and Livine to separate the two directions of propagation in 
the spinfoam sums ^3] can be seen as attempts to separate locally the general relativistic analog of the two 
terms of (j62jl . 



5 Probabilistic interpretation 

The idea developed in |3| is that the variables of the extended configuration space C can be interpreted 
as partial observables, namely quantities that can be measured. A point in C, called an "event", represents 
a simultaneous measurement of partial observables ((a, b) in our case). The classical theory predicts which 
events can be measured: a classical state determines a relation between partial observables and therefore a 
subset of events. The quantum theory assigns a probability amplitude for the measurement of each event. 
According to the general prescription given in [2] , the probability of observing an event in a small region 1Z 
of C if an event has been observed in a small region 1Z' is 

P nn , = \ K KTl>\ 2 ^ 
KfiTl K-R'll' 
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where 

K-jiii 1 = / dadb I dadb' K(a,b,a ,b). (72) 

The propagator K contains then the full physically relevant information about the quantum theory. Alter- 
natively, given a physical state if), the probability for detecting an event around a point (a, b) in a small 
region 1Z of volume V, where ip{a, b) can be taken to constant, is 

Pn = IJ ^y dadb? = ma,b)\ 2 . (73) 

To test the viability of this interpretation in the context of our system, we must relate it to the standard 
interpretation of quantum mechanics. According to the standard interpretation of the wave function i/j(x, t), 
the modulus square of the wave function at fixed time is to be interpreted as the probability density in space 
for the position of the particle. A probability is therefore given by 

dP = \ip(x,t)\ 2 dx. (74) 

Can we recover this expression from (|73|) in a limit in which our system behaves as a nonrelativistic system 
where, say, b is taken as the independent variable (as t) and a as a dynamical variable (as x) varying in 
time? 

For this, we must look for a regime where the fact that the evolution "comes back" in the time b is 
negligible. Consider a region K such that a 2 ,b 2 < M. As in 0, imagine we have a detector active for 
very short "time" b. That is, consider a small detection region Ab <C (Aa) 2 <C 1. Let us then calculate the 
detection probability using (|73|). In the region, r(a, b, a', b') <C 1. Assume also that M S> Ti and we can use 
the saddle point approximation. The propagator takes the form 

, / j /\ 1 . f(a-a') 2 + (b-b') 2 M\ 
K(a, b, a', «) = — sm ^ ^ + -r j (75) 

Using the expression of the parameter r from (|21l) 

yM sin(AVM V(a - aT + ( fe-fe T ) 
K(a, b,a,b) = — ^ a _ d? + {h _ v? ^) 

Since M ^> h, this gives 

K(a,b,af,l/) « ^^y / (a-a') 2 + (6-6') 2 A) (77) 

Performing the integration in the region and using our conditions on the shape of 1Z (which are not 
symmetric in a and b) we obtain 

3 ^ A /A^2 



K nn = f dadb f da'db' K(a,b,a',b') « — — -n- Aa(A6) 2 (78) 



so 



Pr = ]M^ , Mjg H ^ t)fA , (79) 

Kim Aa(Aby 

where we have absorbed a proportionality constant in the normalization of the wave function. The key 
result is that the probability is independent from Ab and proportional to Aa. The first fact allows us to 
use the notion of measurement istantaneous in b, because the "time" b needed by the apparatus to perform 
the measurmenet has no effect on the result. The second implies that \ip(a, b)\ 2 can be interpreted as a 
probability density in a at fixed b. Therefore we have recovered the standard probabilistic interpretation of 
the wave function in quantum mechanics: \ip\ 2 is the probability density in space, at fixed time. 
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6 Conclusion 



We have studied a simple dynamical system that reproduces some key aspects of the background indepen- 
dence of general relativity. We have used this system to illustrate some general features of the structure of 
general covariant dynamical systems. 

According to the point of view we have taken here, the dynamics of a system with n degrees of freedom 
does not describe the evolution in time of n observables. Rather, it describes the correlations between 
n + 1 partial observables. The space of the partial observables is the extended configuration space C, and 
the dynamics is governed by a (vanishing) relativistic hamiltonian H on T*C. In the quantum theory, 
the kinematical Hilbert space /C expresses all potential outcomes of measuments of partial observables. 
Dynamics is a restriction on these states and expresses the existence of correlations among measurements 
of partial observables. Such restriction of the states in /C is given by the Wheeler-deWitt equation Hip = 0. 
Given a physical state, the probability that the system is detected in a small region of C is governed by (|73*|) . 
(A discussion of the precise meaning of probability in this context is in [2].) Here we have shown that in 
our simple system this probability prescription reduces to the conventional one in the appropriate regime. 

This definition of probability yields positive probabilities, unlikely the definition of probability as a 
flux of a current (see Appendix B), often used in quantum cosmology 0. It is reasonable to expect the 
two definitions to agree for states that propagate in only one direction. We also expect it to agree with 
the probabaility computed with histories techniques studied in |Sj and we think that the precise relation 
between these two points of view deserves to be studied. 

The kernel of the projector from the kinematical to the physical state space is the propagator, which 
codes the dynamical content of the theory and can be taken as the basis of the probabilistic interpretation. 
We have studied the propagator of our model in detail. We have shown that in the semiclassical limit it 
has a simple relation with the Hamilton function of the classical theory, but this relation is not a simple 
exponential, as one might have expexted. Instead, the propagator is real. It is the sum of two exponential 
terms complex conjugate to each other, that propagate backward and forward, respectively, along the 
motions. Accordingly, the physical Hilbert space splits between forward and backward propagating states. 
We expect this structure to be the same in quantum general relativity. 
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A Appendix: Approximate Schrodinger picture 

If we want to describe just the local evolution, we can cast the problem in a Schrodinger form by choosing 
an "internal time" variable as a function on the phase-space, before quantizing. The starting point for this 
is to find a function T(a,b,p a ,pt,) that satisfies 

{T,H} = 1, (80) 
We have then dT/dr = 1, here r is the evolution parameter. We chose 

T(b,Pb) = arctan ( — ^ , (81) 
\PbJ 

that satisfies (|80|). Next, we search for a canonnical transformation from the canonical pair b,pb to a 
canonical pair T,Pt- This is given by 



b = ±V2P^smT, (82) 
p b = ± y / 2P^cosT. (83) 

A generating function for this canonical transformation is 
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Indeed, we have 

P = _^ 
dQ 

and 



(85) 



db' 

The Hamiltonian in the canonical variables (a,T,p a , Pt) is then 



P b = ~4- (86) 



H = P T + H eff (87) 

where 

2 



#e// = J(p2 + a 2 )-M 



from which we obtain the Schrodinger equation 

-ih^ = H ef ^ (89) 

It is important to notice that the range of T is contained in the interval [—§,§], therefore this change of 
variables cannto be used around the points b = where the "time" b "comes back". The propagator is, 
up to a phase, the one of the standard harmonic oscillator. 

K(a, T, a, 0) = , 1 exp [ - * - \(a 2 + a' 2 ) cos T - 2a' a] + ~MT} . (90) 
V ^ y/2mh 8 mT \2hsmT L v ' \ h J v ; 

B WKB 

We recall here the interpretation of the semiclassical approximation of a quantum cosmological model in 
terms of the conserved current of the WKB approximation. 

Let us search approximate solutions of the Wheeler-deWitt equation (j.Slj) in the form 

ip{a,b) = A(a,b) e^ s( - a ' b \ (91) 

where A(a, b) is a function that varies slowly, in a sense that we specify in a moment. Inserting this function 
in ()31j) we obtain for the real part 

d 2 A d 2 A A f(dS\ 2 (^\ 2 \ -2(2j + 1) + a 2 + b 2 



wa + m-¥ UrJ + UJ r — i 2 A (92) 



If A(a, b) varies slowly in the sense 



A h 
then S(a, b) must satisfy the Hamilton-Jacobi equation 



8 2 A | 8 2 A rv/r, • I 1 \ I 2 | l2 

WK±WL « -2(2j + l)+a +b (gg) 



1 /5S\ 2 1 /95\ 2 a 2 b 2 fn . 1S n 



For the imaginary part we have 



dAdS dAdS A (d 2 S d 2 S\ . . 

2 ^ +2 ^ +j4 U^ + w =° (95) 
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which can be written as 



d ( A2 dS\ d ( A2 as\ 



da \ da J db \ db 
and interpreted as a continuity equation d a j a + dbj b = for the current j = (j a ,j b ) 

j = A 2 VS. (97) 

where 

- fdS dS\ , . 

VS =U-W < 98 > 

The WKB approximation is a widely used technique in quantum cosmological models. The central question 
is how to extract physical predictions from a WKB solution. The Wheeler-DeWitt equation is typically a 
second-order equation like the Klein-Gordon equation, and the associated current 

J = i(ipVip*-ip*Vtp) (99) 

is conserved: divJ = 0. This current J is intended to provide the interpretation of the WKB solution: the 
flux of J across a surface £ defines the probability that the set of classical trajectories with momentum 
p = VS 1 corresponding to the wave function intersect the surface S. In general, this current produces 
negative probabilities. In particular, we have a full basis of real solutions where J = 0. 



C Path integral representation of the propagator 



Finally, we illustrate a path integral derivation of the propagator. The action of the relativistic system can 
be formulated on an extended configuration space of the coordinates a, b and a lagrangian multiplier N that 
implements the constraint JSJ 



S[a, b, N] 



i / i 

dr I Pa h + p b b - N-{p 2 a +p 2 b + a 2 + b 2 - 2M) 



(100) 



where a dot denotes the derivative with respect to r. The variation of the action with respect to the variables 
a(r),p a (r), b(r),pb(T) and N(r) gives the canonical equations of motion and the constraint equation. The 
action is invariant under reparametrization of the parameter r labeling the motion. This invariance implies 
a relation between the velocities and the momenta, obtained varying the action with respect to p a and 
Pb'.a = Np a and b = Np u . The action can be put in the form 



dr 



'd 2 + b 2 
N 



N(a 2 + b 2 - 2M) 



(101) 



which has the same structure as the ADM action of general relativity. The path integral representation of 
the propagator of the quantized theory has the form 



K(a, b , a, b) 



dN 



Va / Vb e iS[a ' b ' N] 



(102) 



In order to sum over only inequivalent trajectories we fix the gauge in the action choosing N{t) = constant. 
After gauge-fixing and after rescaling the parameter r as f = tN, the propagator can be written 



K(a',b',a,b) = dN Va Vbexp 



a 2 + b 2 - a 2 - b 2 + 2M 
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The range of integration of the lagrangian multiplier will be fixed by the symmetry property of the integrand, 
the explicit expression of the propagator results to be 



K(a , b' , a, b) 



dNVaVbexp 
dN I / £>aexp 



f *!(*>- *)])(/» 



a 2 + b 2 - a 2 - b 2 ) + iMiV 



exp 



v 1 
2 



df-[b 2 -b 2 



AMN 



The expressions in parenthesis have the form of the path integral for the propagator of a free one-dimensional 
harmonic oscillator over the time N in the coordinates a and 6 respectively: 



K(a',b',a,b) = I dNK(a',a,N)K(b',b,N)e lMN 

The propagator of a free harmonic oscillator is a periodic function of the time: 

K(a', a,N + 2vr) = K(a', a, N) 
Consequently we are free to choose the interval [0, 2n] as the range of integration for N 



K 



2- 



dNK(a',a,N)K(b',b,N)e 



iMN 



(103) 



(104) 



(105) 



This representation of the propagator is exactly the expression l|55|). 
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